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INTRODUCTION 


The  determination  of  electron  densities  in  the  middle  atmosphere 

(10-90  km)  is  important  per  se  because  of  the  impact  on  communication 

systems,  and  indirectly  as  an  important  component  in  the  determination  of 

upper  atmosphere  chemistry  .  1  *  2  The  exact  specification  of  ng  in  ionosphere 

models  has  been  a  problem  for  altitudes  below  70  km.  since  most  validated 

3  -3  3  4 

diagnostics  do  not  function  well  with  electron  densities  below  10  cm  .  * 
There  have  been  a  number  of  reported  efforts  which  involve  comparisons  of 
electron  density  indications  in  this  altitude  range  from  more  than  one 
diagnostic  with  data  taken  simultaneously . ^"8  However,  the  results  have 
shown  significant  differences  when  the  data  are  reduced  by  standard  analyses. 
These  analyses  of  diagnostic  interactions  generally  are  based  on  concepts 
developed  in  tandem  with  diagnostic  hardware.  In  many  cases,  there  really 
has  been  no  controlled  verification  of  a  given  ionospheric  diagnostic 
technique,  so  it  is  not  surprising  if  there  are  differences  evident  in 
predicted  composition.  In  fact,  it  is  the  intercomparison  of  reduced  data 
from  diagnostics  that  will  provide  the  most  valuable  test  of  accuracy. ^ 

This  work  will  focus  on  the  definition  of  accuracy  and  resolution  of 
electron  density  by  a  number  of  different  diagnostics.  An  earlier  report"^ 
had  concentrated  on  the  application  of  a  new  analysis  of  electron  collection 
by  blunt  probes  to  data  which  was  taken  simultaneously  with  other  diagnostics. 
One  set  of  tests  used  a  partial  reflection  sounder  and  other  tests  used 
several  rocket  borne  diagnostics.  Reasonably  good  agreement  was  indicated 
in  these  comparisons,  but  there  were  differences.  One  of  these  differences 
was  accounted  for  by  corrections  in  number  densities  derived  from  WSMR 
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partial  reflection  measurements;  these  were  included  in  a  later  reporting. 
The  tests  that  were  evaluated  were  carried  out  in  1972,  1975  and  1977,  while 
the  comparisons  in  electron  density  indications  were  made  in  1979 — well 
after  the  tests.  Upon  summarizing  these  intercomparison  studies,  it  became 
evident  that  only  by  conducting  experiments  with  comparisons  in  mind  before 
the  program  schedules  were  set  could  optimum  intercomparison  of  data  be 
anticipated.  The  1979  US-Canada  Solar  Eclipse  Program  provided  such  a 
campaign  within  which  to  conduct  needed  electron  density  comparisons.  The 
impetus  for  an  intercomparison  effort,  and  indeed  the  work  being  reported 
here,  was  generated  by  R.  Olsen  of  ASL,  WSMR.  The  thrust  of  the  work  to 
be  presented  here  will  be  directed  toward  improving  the  accuracy  of  analysis 
of  data  from  a  number  of  probes  which  were  used  on  that  program:  blunt 
probes;  nose  tip  (supersonic)  Langmuir  probes;  impedance  probes;  and  partial 
reflection  sounders.  The  schedules  of  a  number  of  rocket  firings  during 
the  1979  Solar  Eclipse  were  organized  to  allow  blunt  probe  and  Gerdien  probe 
launches  to  coincide  with  those  of  related  impedance  probe,  tip  probe,  and 
spherical  particle  collector  probe  launches.  This  work  will  be  directed 
to  defining  the  most  exact  analyses  of  data  from  several  of  these  rocket 
borne  probes  such  that  reduction  of  Solar  Eclipse  data  can  be  accomplished 
with  the  most  accurate  procedures.  Efforts  will  focus  on  subsonic  blunt 
probes,  supersonic  tip  Langmuir  probes,  outgas  effects  on  impedance  probes, 
and  shock  or  compression  effects  on  composition  sampled  by  supersonic  probes. 


ELECTRON  DENSITY  DETERMINATION  FROM  CURRENTS 
COLLECTED  BY  SUBSONIC  BLUNT  PROBES 


An  earlier  study  of  blunt  probe  particle  collection  which  resulted  in 

the  prediction  of  electron  density  profiles, ^  exclusively  utilized  a 

technique  of  determining  n  from  conductivity,  O  ,  which  is  proportional 

e 

to  dl/dV,  the  experimentally  determined  current-voltage  slope.  As  noted  in 

12 

those  works,  however,  it  is  possible  to  apply  a  new  analysis  of  the 

particle  collection  to  derive  electron  densities  directly  from  currents 

collected  at  specific  voltages,  for  a  given  probe  shape  at  a  given  altitude. 

That  technique  will  be  applied  here  to  data  that  was  gathered  with  two 

different  types  of  probes:  subsonic  blunt  probes  and  supersonic  nose  tip 

probes  with  approximately  conical  shape.  As  in  the  earlier  work,  the 

indications  of  electron  density  that  are  derived  will  be  compared  with  those 

deduced  from  conductivity  and  with  those  from  other  diagnostics. 

For  the  subsonic  blunt  probe,  several  new  pieces  of  information  can  he 

added  to  the  present  understanding.  First,  it  will  be  noted  that  the 

4  5 

standard  procedure  adopted  by  the  Illinois  group  ’  which  gathers  electron 
saturation  currents  by  biasing  the  tip  electrode  at  a  constant  2.7  V,  is 
to  "calibrate"  the  probe  on  a  given  launch  by  matching  the  current  to  an 
electron  density  at  one  altitude.  The  variation  of  electron  density  with 
altitude  is  thus  implied  by  the  variation  of  electron  saturation  current. 

As  noted  in  an  earlier  work,  this  may  have  some  validity  at  higher  ( > 70  km) 
altitudes,  but  the  procedure  leads  to  questionable  implications  of  np  at 
lower  altitudes  because  of  inter-particle  collision  effects  on  collection, 
along  with  shock  wave  and  electron  reattachment  effects.  These  effects 
have  not  been  specifically  calculated  for  any  such  probe. 
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Again,  accepting  the  fact  that  tip  probe  data  has  traditionally  been 
presented  as  I  (current  at  2.7  V)  and  blunt  probe  data  has  traditionally 
been  presented  as  0_  (dl/dV),  it  is  of  interest  to  evaluate  and  to  graph 
subsonic  blunt  probe  electron  saturation  currents  to  compare  with  simulta¬ 
neous  supersonic  tip  probe  currents.  This  comparison  has  been  able  to  be 
carried  out  with  data  that  was  gathered  by  Hale  (Penn  State  University) 
simultaneously  with  measurements  taken  by  the  Illinois  group  during  the 
Winter  Anomaly  campaign".  Data  was  available  for  January  31,  1972  and 
December  5,  1972.  For  the  Super  Areas  configuration,  where  a  (probe  radius) 
3.65  cm,  the  currents  collected  by  the  blunt  probe  were  evaluated  at  V  =2  V, 
and  were  graphed  in  the  Illinois  fashion,  i.e.,  "calibrated"  by  relating  the 
current  to  electron  density  at  a  given  altitude.  On  Figure  1,  data  from 
January  31,  1972  are  presented;  Illinois  data  are  as  published,  where  the  tip 
probe  data  are  calibrated  by  the  use  of  Faraday  rotation  and/or  differential 
absorption  rocket-borne  techniques.  For  reference,  n  Indications  determine 
from  the  a_  reduction  are  shown.  Since  currents  from  the  blunt  probe 
at  2  V  were  not  available  at  altitudes  higher  than  65  km,  the  current  was 
"calibrated"  to  the  Illinois  current  at  that  point.  Based  only  on  the 
comparison  of  I  (blunt)  and  ng  (o_) ,  there  are  no  great  disparities  in  the 
relative  variation  with  altitude.  This  same  procedure  was  carried  cut  for 
the  data  available  from  December  5,  1972;  this  is  presented  in  Figure  2. 

In  this  instance,  it  was  possible  to  "calibrate"  the  blunt  probe  current  at 
the  highest  altitude  (74  km)  with  density  from  more  accurate  rocket  instru¬ 
ments.  Generally  speaking,  the  agreements  with  Illinois  I  is  good  from 
74  km  down  to  60  km  where  it  ends.  It  is  worth  noting  that  the  unique 
reversal  in  ng  at  64  km  on  this  shot  is  mirrored  on  both  blunt  and  tip 
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density  variations  for  January  31,  19/2,  1230  tST,  Wallops  Island: 


Figure  ?.  Electron  density  variations  for  December 
blunt  probe  Ip  and  calculated  densities 


probes.  From  the  above,  it  on;,  oe  concluded  that  there  is  little  apparent 

difference  oi  I  variation  provided  by  supersonic  tip  probes  and  subsonic 

blunt  probes  at  altitudes  above  60  km. 

With  the  electron  saturation  current  at  a  given  voltage  categorized 

for  the  blunt  probe,  it  is  now  of  interest  to  determine  the  electron  density 

based  on  I  at  V  data.  For  the  blunt  probe,  this  is  straightforward,  as  the 
e  P 

relevant  ecuations  had  already  been  established ^  For  the  two  days 

noted  above,  the  appropriate  equations  and  reduced  electron  densities,  n  ,  arr 

also  presented  in  Figures  1  and  2.  It  will  be  emphasized  at  this  point,  t ha' 

the  n  indications  are  carried  out  without  anv  correction  or  calibration: 
e  — 

there  is  no  lateral  shift  to  effect  agreement.  The  n&  values  in  these 

figures  exhibit  the  same  order  of  magnitude  as  those  from  other  diagnostic':: 

the  relative  variations  are  also  similar.  One  very  interesting  aspect  or 

these  comparisons  can  be  noted  here:  for  December  5,  1972  the  I  indications 

e 

are  less  than  0  indications  of  n  ,  while  for  January  31,  1972  the-  I  indi- 

e  e 

cations  are  greater  than  C_  indication  of  ng .  One  difference  between  T 

indications  and  0  indications  is  that  I  can  be  perturbed  bv  local  E  fields 

e 

or  any  variation  of  vehicle  potential  from  reference  potential.  Accordingly. 

such  potential  or  field  differences  can  be  evaluated  on  that  basis. 

These  same  procedures  can  be  applied  to  data  taken  with  a  Loki-Dart 

sized  blunt  probe,  where  a  (p-obe  radius)  =2.08  cm,  at  White  Sands  Missile 

Range  (WSMR)  on  two  different  days.  As  was  done  earlier,  the  blunt  probe 

indications  of  n^  can  be  compared  with  partial  reflection  data.  Partial 

reflection  and  n  derived  from  C_  can  be  considered  reference  data  in 
e 

Figures  3  and  4.  Electron  current  collected  with  V  =  2  V  is  graphed  in 

P 

these  figures  with  a  "calibration"  point  chosen  at  highest  altitude. 


11 


(  BLUNT  ) 


Figure  3.  Electron  density  variations  for  October  1975,  1 6 >  'DT,  WSMR:  blunt  probe 
and  calculated  densities  from  I„. 


’  variations  for  September  ?9,  19/7,  1345  4ST,  WSMR:  blunt 
culated  densities  from  I  . 


Electron  densities  determined  by  exact  formula,  without  lateral  shaft,  are 

also  presented  in  the  figures.  As  for  the  data  presented  from  Wallops 

Island  during  1972,  there  are  similar  orders  of  magnitude  and  variations 

with  altitude  from  all  diagnostic  techniques  and  procedures.  Also,  as 

above,  the  1977  day  shot*  n  from  I  below  n  from  0_,  while  the  1975  day 

e  e  e 

shows  n  from  I  almost  exactlv  the  same  as  n  from  Q_  for  48-60  km.  This, 
e  e  e 

again,  could  be  caused  by  potential  shifts  or  local  field  effects. 

ELECTRON  DENSITIES  DETERMINED  FROM  SUPERSONIC  TIP 
DATA — COMPRESSION  AND  ELECTRON  ATTACHMENT  EFFECTS 

The  above  comparisons  were  based  on  published  information  and  subsonic 
blunt  probe  data  which  had  been  shown  in  an  earlier  work  to  be  somewhat 
comparable.  However,  one  of  the  most  important  tests  of  the  analysis  of 
data  from  electron  collecting  probes  must  come  from  different  geometry 
probes,  flown  by  different  investigators.  Evaluation  of  such  data  to  deter¬ 
mine  n  ,  and  the  comparison  of  n  with  that  already  presented  would  be 
e  e 

important  in  establishing  the  universality  of  any  technique.  Fortunately, 
raw  data  gathered  by  the  Illinois  group  using  supersonic  nose  tip  Langmuir 
probes  during  the  1972  days  at  Wallops  Island  were  available.  This  data  was 
retrieved  and  provided  by  L.  Smith.  Also,  the  detailed  geometry  oc  the  tip 
collector  was  specified  as  being  an  ogive  with  height  (base  to  tip)  of  1.5  in 
base  diameter  of  1.082  in.  and  an  arc  radius  of  curvature  of  3.209  in.  While 
calculations  can  be  carried  out  with  reasonable  accuracy  presuming  a  conical 
shape  instead  of  an  ogival  shape,  in  this  work  a  further  approximation  was 
made:  the  ogive  was  represented  by  a  hemisphere  of  equivalent  area.  A  hemi¬ 

sphere  of  radius  1.42  in.  would  have  equal  area  with  the  Illinois  ogive. 
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Following  established  procedures,  the  electron  flux  collected  by  a  hemi¬ 
sphere  can  be  written 


n  v 
e“>  DB 


2Trr; 


1 

4  nXC 


2ir(r 


A  )‘ 
en 


where  n  is  the  undisturbed  electron  density,  is  the  saturation  velocity 

e°°  DB 

at  r_ ,  r  =  5  r  (Ew/p)  is  the  radius  for  saturation  flux,  n,  is  the  electron 

D  D  p  ), 

density  at  one  mean  free  oath,  A  ,  from  the  collector  surface,  c  is  the 

en 

electron  thermal  velocitv,  and  r  is  the  radius  of  the  hemispherical 

P 

collector  (probe)  surface.  Incorporating  this  into  a  formulation  for  curren' 
collected  by  the  probe, 


I 

e 


e[-^  n^c](27rr2) 


results  in 


the  following  formula  for  electron  density  for 


r 

P 


V  =  2.1  V: 
P 


1.81  cm. 


10  2 
n  =  (2.72  x  10  )I  (A)p(mm  Hg) (1  +  A  /r  ) 
e°°  e  r  a  en  p 


When  appropriate  data  were  substituted  into  the  above  formula,  the  values 
of  electron  densities  were  calculated  in  a  straightforward  manner.  These 
electron  densities  are  indicated  in  Figures  5  and  6  (bv  the  notation  "N  ‘'rot 
Supersonic  Tip  Probe");  generally  these  values  are  higher  at  higher  altitude 
and  lower  than  other  diagnostics  at  lower  altitude.  At  altitude  above  60  km, 
the  rate  of  change  of  n^  is  in  reasonable  agreement  with  that  derived  from 
other  diagnostics.  Before  proceeding  with  further  calculations,  it  should 
be  noted  that  the  analytical  formula  above  is  sensitive  to  the  model  chosen 
for  A^,  the  collection  area  at  A  from  the  probe  surface.  However,  using 
a  consistent  approach,  an  alternative  approximation  with  cylindrical  shape 
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pave  the  same  numerical  results.  Accepting  the  above  form  as  reasonable, 
further  analyses  will  proceed. 

There  are  two  effects  that  must  receive  attention  in  the  reduction  of 

supersonic  collection  data:  (1)  the  compression  effects  of  shock  waves 

about  the  collecting  shape,  and  (2)  the  alteration  of  plasma  chemistry  by 

the  shock  wave  event.  Both  of  these  effects  are  dependent  upon  the  vehicle 

trajectory  (Mach  number).  It  will  be  reiterated  that  an  earlier  report 

identified  a  number  of  references  which  established  the  magnitude  and  extent 

oc  shock-type  compressions  un  to  100  km.  Further,  it  can  be  roted  that  the 

electron  collection  processes  generate  particle  flux  to  the  surface  at 

radius,  r  ,  which  is  far  from  the  surface;  r  will  be  well  outside  any 
B  B 

shock  or  compression  region.  However,  kinetic  flux  to  the  surface  at  1 

will  be  proportionally  altered  by  the  shock  or  compression  effects.  A 

simple,  linear  approach  allows  the  compression  effects  to  be  compensated 

for  in  a  seauential  manner.  In  order  to  evaluate  the  magnitude  of  these  two 

effects  or  electron  density,  approximate  calculations  were  carried  out. 

Firs-,  based  on  trajectory  information  supplied  by  Illinois,  Mach  numbers 

were  evaluated  as  a  function  of  altitude.  Or  January  31,  1972.  flight  Mac'. 

number  was  about  3.2  for  50-80  km;  on  December  5,  1972,  flight  Mach  number 

was  about  4.9  for  50-80  km.  For  the  ogive  shape  used,  the  half-cone  angle 
1.  3 

of  30°  is  indicated  to  result  in  a  shock  wave  angle  of  about  35°.  By 
standard  procedures,  the  density  ratio  across  the  shock  wave  is  P  /?.,  -  4 
and  T2/T1  ~  2.7.  Directly,  this  means  that  the  density  of  particles  with¬ 
in  the  shock  cone  would  be  raised  by  a  factor  of  4;  to  recover  densities 
in  the  ambient  medium,  it  is  appropriate  to  divide  the  indicated  electron 
densities  by  this  factor  of  4.  The  primary  values  of  electron  density  noted 
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above  were  adjusted  first  for  shock  compression  effects,  and  these  are 

noted  in  Figures  5  and  6  by  arrows  indicating  the  adjustment.  It  can  be 

seen  that  this  calculation  procedure  results  in  values  that  are  in  very  close 

agreement  with  those  of  other  diagnostics  for  altitudes  between  80  and  60  km. 

Below  60  km,  the  electron  densities  reduced  by  including  shock  compression 

effects  are  substantially  lower  than  those  derived  from  subsonic  blunt  probes. 

A  second  major  correction  effect  to  the  supersonic  nose-tip  Langmuir 

probe  data  must  be  considered  at  lower  altitudes.  The  presence  of  a  strong 

compression  region  in  front  of  the  collector  can  result  in  substantial 

electron  loss  through  enhanced  reattachment  of  the  electrons.  Specif ical lv . 

1A  i  s 

electron  loss  mechanisms  were  reviewed  by  both  Mitchell  and  Heaps,  and 
the  dominance  of  loss  by  attachment  was  identified.  The  simplest  model  that 
can  describe  attachment  is  written 


dn 
_ e 

dt 


^  -  8  n 


^here  8  is  the  attachment  rate  coefficient  derived  from  the  two  reaction* 


and 
to  be 


N'2  +  +  e  -*■  0,  + 

°2  +  °2  +  €"  °2  +  °n 


.  -7.9  .300.  -(boon) 

l.A  x  10  (-=-)e 


|0  T  +  10 


-31, 


0  J  'V 


In  the  above,  [  ]  refers  to  the  number  densitv  of  the 


thesis.  Integration  of  the  equation  for  n?  gives 


-Pt 

n  =  n  e 
e  eo 


in  the  ;  mo 
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spec  i  e s 


where  n  is  the  electron  density  behind  the  shock  wave  after  chemical 
e  - 

reaction  involving  incoming  flow  has  taken  place  while  ngo  would  be  the 
density  behind  the  shock  waves  presuming  no  chemical  reaction  (reattachment). 
The  characteristic  time  available  for  reattachment  can  be  evaluated  by  con¬ 
sidering  the  average  speed  of  particles  behind  the  shock  and  the  average 
distance  traveled  before  collection; 

4  -1  -4 

t  ~  d/V  -  1  cm/10  cm  sec  =  1C  sec  . 

Accordingly,  the  effects  of  reattachment  are  included  by  using  the  shock 

adjusted  values  of  n  noted  above,  and  aDprooriate  8,  t  values  to  calculate 

e 

n  .  It  should  also  be  noted  that  the  value  of  8  that  was  used  is  different 
e 

from  that  appropriate  for  ambient  conditions,  since  the  shock  compression 

and  heating  alter  the  density  and  temperature  values  and  also  B.  Further, 

t^ere  is  some  conjecture  regarding  'he  correct  attachment  rates  in  the  light 

of  available  information.1^  Specifically,  Heaps"1^  notes  that  attachment 

rates  about  two  orders  of  magnitude  higher  than  standard  rates  are  necessary 

to  allow  his  computational  model  to  he  in  reasonable  agreement  with  changes 

noted  during  an  eclinse.  Tn  ‘‘act,  this  same  type  of  behavior  is  consistent 

with  supersonic  tip  probe  data:  if  an  attachment  rate  is  chosen  which  is 

a  factor  o f  li>  higtier  than  standard,  there  is  good  agreement  between  subsonic 

■  '  unt  pro'll'  n  and  supersonic  '  ip  probe  n  at  all  altitudes  below  60  km. 

>•  e 

The  appropriate  data  points  are  shown  in  Figures  5  and  6,  designated  with 
a  second  arrow  shift  to  higher  values  to  correct  for  electron  attachment. 


20 


Based  on  all  of  the  above  considerations.  It  Is  of  Interest  to  compare 
the  most  accurate  indications  of  ng  from  blunt  probe  o_,  with  ng  from  Illinois 
nose-tip  Ie  data,  and  Illinois  Faraday  rotation  and  differential  absorption 
points  above  70  km.  These  are  shown  in  Figures  7  and  8.  While  Figure  8  shows 
good  agreement  at  all  altitudes,  this  was  an  abnormal  day.  Figure  7  shows  tip 
probe  indications  of  ng  from  Ig  below  60  km  to  be  lower  than  blunt  probe  ne 
for  a_;  it  is  felt  that  the  blunt  probe  technique  and  analysis  is  most 
accurate  here.  However,  above  60  km,  the  blunt  probe  indications  of  ng  from 
a_  are  significantly  lower  than  other  diagnostics.  Blunt  probe  methods  above 
65  km  are  approximate,  and  it  is  possible  that  these  lower  indicated  densities 
are  due  to  the  approximate  method  of  analysis  used.  These  will  be  further 
explored  in  the  future. 

OUTGAS  PERTURBATIONS  OF  MEDIUM  SURROUNDING  ROCKET  BORNE  PROBES 

The  primary  concern  in  this  work  is  with  the  intercomparison  of  electron 
density  indicative  diagnostics  in  the  middle  atmosphere.  Particle  collecting 
subsonic  blunt  probes  and  supersonic  tip  probes  have  been  given  detailed  con¬ 
sideration,  and  the  results  of  analyzing  their  data  have  been  compared  with 
other  diagnostics.  Generally  speaking,  those  diagnostics  are  not  sensitive 
to  the  presence  of  outgas  material,  since  they  function  as  collisional 
particle  collectors  and  the  mobility  of  electrons  in  outgas  species  would 
not  generally  inhibit  collection.  However,  there  are  diagnostics  that  could 
be  sensitive  to  the  presence  of  outgas  species,  and  so  it  is  important  to 
categorize  the  amount  of  outgas  and  the  spatial  distribution  of  outgas 
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Figure  7.  Electron  density  variations  for  January  31,  19/?,  I?30  K ST,  Wallops  Island: 

from  blunt  probe,  supersonic  tip  probe  and  local  wave  interaction  methods. 
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materials.  In  particular,  one  of  the  most  important  diagnostics  for 
intercomparison  purposes  uses  RF  probe  techniques.^  The  theory  of  these 
diagnostics  has  shown  a  dependence  on  electron  collision  frequency,  and  so 
consideration  of  outgas  effects  will  be  directed  to  that  application.  The 
work  to  be  presented  in  this  section  will  be  a  condensation  of  work  that  was 
informally  reported  earlier.^ 

Of  particular  interest  here  is  the  renorted  observation  bv  Kendall 

that,  under  the  conditions  of  explosive  decompression  of  the  atmosphere 

surrounding  a  rocket-borne  ionosphere  nrobe,  the  rate  of  outgassine  from  the 

probe  appears  to  be  orders  of  magnitude  larger  than  previously  anticipated. 
i  8 

Standard  estimates  of  outgas  effects  use  rates  that  are  more  representa¬ 
tive  of  conditions  after  pumping  periods  of  hours  or  hundreds  of  hours. 
Accordingly,  the  comoosition  of  the  gases  near  the  surface  of  rocket-launched 
probes  with  a  flight  period  of  5  to  10  minutes  can  he  expected  to  be 
significantly  different  from  earlier  estimates.  Since  the  gas  desorbed  by 
probe  components  is  primarily  water  vapor,  a  natural  component  of  the  iono¬ 
sphere  whose  density  and  chemical  role  is  not  well  understood,  any  such 
species  contamination  is  of  considerable  interest. 

An  attempt  will  be  made  here  to  estimate  order-of-magni  tude  values  >  •' 
the  water  vaoor  density  profiles  near  a  desorbing  probe  surface  so  as  to 
establish  the  extent  of  any  distortion  of  ambient  gas  composition  by  * 
probe  itself.  The  variation  of  concentration  as  a  function  of  distance  from 
a  plane,  outgassing  surface  will  be  considered. 

The  problem  of  determining  the  distribution  of  outgassin g  molecules 
near  a  solid-gas  boundary  involves  consideration  of  the  mass  transfer  at  the 
boundary  and  the  subsequent  species  diffusion;  it  possesses  some  similarity 
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i  y  .  ' 

to  the  problem  of  evaporation  (liquid-gas)  and  ablation  (so  1 i d- 1 i q ui d- 

gas)  in  high  temperature  flows.  The  analysis  of  those  two  problems  involves 
a  coupling  between  mass  and  energy  transfer,  along  with  the  required  momentum 
balance  in  the  situation  with  flow.  The  constraint  normally  applied  at  the 
boundary  with  the  gas  is  the  specification  of  equilibrium  between  the 
injected  gas  species  (A)  and  the  original  gases  present,  i.e.,  the  partial 
pressure  of  suedes  A  is  the  vapor  pressure  corresponding  to  the  ndxtv.r' 
pressure  and  temperature.  Such  a  boundary  condition  is  net  meaningful  ir 
the  outgas  situation  where  the  rate  of  molecules  being  introduced  at  the 
solid-gas  surface  is  controlled  by  the  mechanism  of  pressure  imbalance 
between  the  surrounding  gases  and  the  adsorbed  species,  and  is  significant!'.' 
influenced  by  the  local  surface  physics. 

A  review  of  current  understanding  of  the  mechanisms  involved  in  out- 

21 

gassing  is  presented  by  Santeler  et  al .  Briefly,  a  gas  mixture  in  contact 

with  a  solid  surface  will,  after  a  long  period  of  time,  reach  a  stite  of 

equilibrium  where  the  amount  of  gas  striking  and  adhering  to  the  surface 

will  equal  the  amount  of  gas  leaving  the  surface.  Absorption  (surface'1  and 

absorption  (volume)  equilibrium  will  be  reached.  One  possible  model  for  t'"o 

determination  of  surface  partial  pressure,  while  clearlv  approximate,  mav 

serve  to  indicate  the  surface-gas  interaction  in  the  outgas  process. 

In  explosive  decompression  it  may  be  assumed  that  loss  of  surface 

particles  from  adsorption  will  occur  first.  Accordingly,  q,  (0),  the  efflux 

2 

of  molecules  (torr-liter/cm  -sec)  from  the  surface  at  t  =  C  can  be  written 


S  _ 

c,  (0)  =  --  =  cxP  (3.638)  /T/MW 
*L  T_  o 
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where  =  quantity  (number  of  molecules)  of  gas  residing  in  a  monolayer 
on  the  surface 

x  =  residence  time  of  molecules  in  a  monolaver  on  the  surface 
R 

(10  -  100  layers),  sec 

a  =  sticking  fraction  for  molecules  impinging  on  the  surface 
P  =  partial  pressure  of  vapor  in  equilibrium  at  t  =  0,  torr 

O' 

T  =  temperature,  °K 

M  =  molecular  weight  of  the  gas,  gr./jrm-m ole 


The  above  equation  relates  atomic  surface  properties  (F^,  "  ,  ex'  v"  th  ear- 
properties  (P  ,  T,  MW).  It  is  presumed  that  these  (material)  surface 
properties  can  be  determined.  Further,  after  the  initiation  of  desorption, 

if  it  is  assumed  that  the  outgassed  molecules  will  be  lost  at  such  a  rate 
so  as  to  establish  a  quasi-steady  equilibrium  with  the  wall,  then 

S  j  T 

qT  (t )  =  e  -R  =  P(3.638)  /tTmT  . 

This  provides  a  possible  scheme  for  determining  the  effective  near-surf, ?• r 
partial  pressure  (?)  of  the  outgassed  species  as  a  function  of  time.  V."-  -  1  •  ■ 
the  accuracy  of  the  above  procedure  clearly  rests  upon  a  precise  knowledge 
of  surface  properties,  it  does  serve  to  emphasize  the  relative  indeoenderce 


of  the  surface  boundary  condition  from  either  diffusion  and/or  flow  ec.' 


Models  which  preserve  the  dominance  of  such  effects  will  now  be  examined. 


TRANSIENT  DIFFUSION  INTO  A  SEMI-INFINITE  MEDIUM 


The  physical  problem  being  considered  -’nvolves  a  solid  with  adsorbec 
gas  being  subjected  to  a  surrounding  gas  pressure  which  will  decreases 
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from  sea  level  atmosphere  to  effective  vacuum  over  approximately  a  three- 


minute  time  period.  The  dynamic  development  of  water  vapor  concentration 

profiles  will  be  determined  in  a  model  where  a  plane  solid  is  subjected  to 

a  vacuum  AP  increment  at  t  =0.  The  analysis  of  such  a  system  is  similar 
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to  that  for  evaporation  of  a  liquid  in  a  tube  of  infinite  length  with  1-D 
diffusion.  The  basic  solution  is  carried  out  for  the  case  where  the  gas  is 
maintained  at  constant  pressure  and  temperature  after  initiation  (t  =  9). 

The  solution  is  carried  out  in  terms  of  dimensionless  variables,  as  there  is 
no  characteristic  length  to  be  defined.  With  zero  outgas  concentration  at 
t  =  0  and  z  (normal  distance)  =  °°,  the  variation  of  outgassed  species,  is 
given  by 


-  =  1  -  erf(Z  -  <t>) 
1  +  erf  (J> 


where  X 


XA/XA0  ’  XAO 


Vz 


0) 


Z  -  z(4Dab  t)-1/2 

Y  _  ^A  r  =  _^A  .mole. 

XA  C  ’  A  MW  ^  V 
A  cm 


z  =  distance  from  solid-gas  interface 


e  =  i/c  Stpo 


AB 


N 


3 

C  =  mixture  concentration  [C  =  p/RT  =  d/MW  (moles/cm  )] 
t  =  time  (seconds) 

D. „  =  coefficient  of  mutual  diffusion 
AB 

2 

N  =  flux  of  particles  at  z  =  0  (gm-moles/cm  -sec) 

Ao 

Thus,  with  a  formula  for  D  ,  specification  of  p,  T,  z,  t,  N»  (outgas  flux'1 

AB  ° 

will  allow  evaluation  of  accordingly,  the  profile  of  X  /XA  can  be 

determined.  However,  since  no  boundary  condition  or  composition  has  been 

specified,  both  X,(z,t)  and  XA  (t)  are  unknown. 

A  nO 
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It  should  be  noted  that  in  the  form  presented  above,  the  solution  is 


exact  only  for  small  values  of  X^  i.e.,  where  outgas  species  A  is  diffusing 
through  a  background  species,  B.  In  order  to  place  the  solution  in 
perspective,  one  calculation  of  interest  involves  the  determination  of  the 
distance  from  the  wall  within  which  90  percent  of  the  outgas  mass  is 
contained  i.e.,  where  X^/Xaq  =  0.1;  for  convenience,  we  call  this  distance 
z(0.1).  As  a  simplifying  assumption,  a  linear  variation  of  composition  is 
utilized  in  some  integrations  to  arrive  at  a  closed  form  expression, 

z(0.1)  =  2.3  /D~~  t1/2  +  2  —2-  t 

AB  C 


Following  Bird,  Steward,  and  Lightfoot,^ 


rcm2) 


=  2*13/2  r_l_  +  1  1 1/2 


,3/1 


2m, 


p(dA  +  dB/2) 


2  2  —8  7 

so  at  50  km,  D„ „  =  4.1  x  10  cm  /sec,  and  C(50  km)  =  3.5  x  10  gm-mole/cm. 
Ad 

Using  the  outgas  rate  indicated  by  Kendall^7 (worst  case,  anodized  aluminum) 
—8  2 

N  =  1.7  x  10  gm-moles/cm  -sec,  then  for 
Ao 


-2 

t  =  10  sec  ,  z(0.1)  =  4.6  cm 

t  =  1  sec  ,  z(0.1)  =  48  cm 


Diffusion  of  outgas  species  will  proceed  to  a  steady  state  distribution  in 
a  time  period  short  compared  with  the  characteristic  time  of  the  local 
probing  event  (y  1  sec).  Alternatively,  it  can  be  stated  that  the  diffusion 
length  scale  will  considerably  exceed  the  scale  of  the  probe  geometry. 
Accordingly,  in  the  present  case  the  diffusion  can  be  treated  as  steady 
state  and  so  will  be  controlled  by  probe  and  flow  generated  dimensions  and 
boundary  conditions. 
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STEADY  STATE  DIFFUSION  THROUGH  A  BOUNDARY  LAYER 


The  model  chosen  is  normally  used  to  represent  diffusion  through  a  tube 
from  which  particles  are  carried  away  (X^  =  0)  at  one  end  (z  =  z )  by  a  gas 
flow  perpendicular  to  the  direction  of  diffusion.  Such  a  model  can  be  used 
to  represent  the  diffusion  of  outgassed  vapor  in  a  case  with  flow  impinge¬ 
ment,  if  an  appropriate  length  scale  can  be  selected  to  represent  the  distance 
from  the  boundary  within  which  the  flow  will  not  disturb  the  diffusion.  Such 
a  dimension  can  be  defined  and  related  to  a  boundary  layer  thickness. 


19 

As  outlined  in  the  literature  the  diffusive  flux  of  A  particles 
through  medium  B  can  be  expressed  by 


Vz) 


CD,,,  dx* 
AB  _ A 

L  -  X,  dZ 
A 


X  =  — 
*A  C  ’ 


C  =  — - 
A  MW, 


with  the  boundary  conditions 


N 

outgas 


cm  -sec 


0 


X.  =  0  at  z  =  z  , 

A  2 

This  equation  can  be  integrated  to  predict  the  concentration  ratio  at  the 
surface  as 

X  =  i  -  e~*NAo/DAB  *  C}z2 
A 

o 

Determination  of  outgas  species  distribution  in  the  case  with  flow 
impingement  presents  a  complex  problem.  The  velocity  flow  field  can  be 
separated  into  two  regions:  (1)  an  inviscid  (external)  flow  region  and 
(2)  flow  within  a  viscous  boundary  layer.  The  thickness  of  the  boundarv 
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layer  presents  a  natural  length  with  which  to  scale  the  flow  problem. 

14 

Specifically,  with  laminar  flow  along  a  flat  plate,  and  the  outgas  <1 

species  A  at  the  plate  with  concentration  t'A  ,  the  A  species  will  diffus« 

no 

and  have  non-zero  concentrat ion  within  a  diffusion  layer  of  thickness,  ' 
approximately  specified  as 

a  =  *J>  *  (Scr1/H  .  r  1/3 


The  absence  of  chemical  reactions  is  presumed.  ''value*  ;  nv  fee  a  Is  ••  •r.i-*- 

eter  for  a  specific  altitude,  A  ( Y.  km'  =  l.?1,  so  t'-*-  v '  •  •  vtu1'.  v  • 

thickness  provides  a  reasonable  representat  ion  of  the  region  w.  » *>  < «  v- 

outgas sed  molecules  will  he  distributed  along  a  c 1  a t  plate.  V-  i  ’  • • 

stagnation  region  of  a  flow  normal  to  a  flat  plate  is  of  some  interest  here, 

these  results  for  flow  along  a  flat  plate  also  provide  a  reievar*  'eng:1 

scale  for  stagnation  boundary  flow. 

Fundamentally,  the  flow  field  near  the  approximatt  lv  plane  boundary 

will  alwavs  be  subsonic  in  nature,  since  even  with  supersonic  flow,  a  normal 

shock  wave  will  convert  the  incoming  supersonic  flow  to  su'-'s  nic  flow,  wit 

a  density  higher  than  ambient  behind  the  shock.  The  incomress ! u le  stacna- 

o  7 

tion  region  is  considered  by  Schlichting*""  where  the  thickness  o r  the  axi- 
symmetric  boundary  layer  is  given  as 


6  =  2.8  /v/ 2a  where  a  =  du/dr  ,  v  =  u/d 

For  a  flat-nosed,  right-circular  cylinder 


a 


4U 

D 


with  U  =  free  stream  velocity  and 
D  =  diameter  of  cylinder 
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ami  d : recti v,  with  L'  -  100  m/s ec,  0  =  0.1  m, 

6(50  km)  =  0.65  cm 

6(70  km)  =  2.0  cm 

6(90  km)  =  11.0  cm 

With  an  appeal  to  the  "similarity"  of  flow  within  the  boundary  lavor 
of  tangent  and  norma!  flat  plate  flows,  the  approximation  that 
t  inr  o!  '’it  gas  molecules  will  be  effectively  zero  at  the  outer  ere.  ■  *  • 

diffusive  and  viscous  boundary  layer  thickness  can  1  ■  madt  ,  in  ’  a*>  or  :>  r 
magnitude  for  the  surface  concentration  of  the  outnas  ;  e<  ui”  arr 

at.  Presuming  steady  state  diffusion  with  no  mass  huilo-uo  t  - rough  a 
boundary  laver  thickness,  the  relationship  for  concent  rat  ion  of  out  <-: 
at  t  he  sur  face  is: 

=  1  _  €'(NA0/cdAB  •  -■  D 
C 

_g  2  o 

Wita  Na  =  1*7  x  10  go-tnoles/cm  -sec,  6  *  £•  ,  then  CA  /C  10”  .  Ar  50 
°  D  no 

the  indicated  density  of  water  vapor  at  the  surface  is  nHo0  =  2  x 10  /cirh 
one  to  two  orders  of  magnitude  higher  than  the  water  vapor  densitv  in  the 
ambient,  undisturbed  atmosphere. 

Based  on  the  above  calculations,  consideration  will  now  be  given  to  t 
influence  of  outgas  water  vapor  on  electron  density  diagnostics  wit1'  the 
impedance  probe. 

INFLUENCE  OF  OUTGAS  SPECIES  ON  THE  REDUCTION  OF  ELECTRON 
DENSITIES  FROM  IMPEDANCE  PROBE  MEASUREMENTS 

The  use  of  RF  probing  techniques  for  the  determination  of  electron 
densities  has  been  accepted  over  a  period  of  years.7  One  early  work 
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satisfactorily  compared  the  general  performance  behavior  of  three  different 

techniques:  (1)  the  standing  wave  impedance  probe  or  impedance  probe; 

(2)  the  plasma  frequency  probe  or  plasma  probe;  and  (3)  a  resonance 

rectification  probe.  'Hie  analysis  of  the  interactions  occurring  with  each 

of  these  types  of  diagnostics  was  reviewed  in  some  detail  by  Despain.^3  An 

attempt  at  quantitative  comparison  of  n^  from  impedance  probes  and  the  other 

24 

diagnostics  was  reported  by  VI wick  in  relation  to  data  gathered  curing 

a  solar  particle  event  (PCA  69).  However,  these  n^  indications  were  not 

absolute  and  the  altitudes  of  best  accuracy  were  from  80  km  to  125  km.  An 

attempt  at  intercomparison  of  an  RF  probe,  the  plasma  frequency  tyee,  with 

25 

a  DC  Langmuir  probe  was  carried  out  by  Kotadia  and  Kist  who  reported 
reasonable  agreement  of  measurements  taken  from  95  to  150  km. 

Basic  analyses  of  RF  interactions  with  plasmas  have  been  presented  by 

a  number  of  authors,  and  in  some  cases  experiments  have  been  reported. 

26 

Specifically,  Balmain  presented  a  basic  theory  and  some  confirming 

laboratory  experiments  with  a  cylindrical  antenna.  Sheath  effects  and 

collisionswere  more  accurately  modeled  in  a  later  work,"'7  but  in  general, 

most  interest  was  focused  toward  plasma  frequency  diagnostics.  A  spherical 

impedance  probe  was  evaluated  in  a  laboratory  experiment  by  Tarstrup  and 

Heikkila,^®  who  examined  fluid  effects  and  effects  of  electron-neutral  colli 

sion  frequency.  Again,  conditions  near  plasma  frequency  were  most  important 
29 

A  later  work  examined  transient  sheath  effects.  A  reasonably  comprehen¬ 
sive  analysis  of  the  impedance  of  a  spherical  probe  in  a  warm  plasma,  with 
collision  effects,  was  presented  by  Aso.^®  In  almost  all  of  the  works 
noted  above,  the  emphasis  was  directed  toward  exploring  behavior  near 
plasma  frequency. 
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A  somewhat  different  emphasis  is  involved  in  the  impedance  probe  ap 

tion  as  outlined  by  Despain.^3  The  impedance  probe  technique  general Iv 

functions  with  to,  the  operating  or  driving  frequency,  being  lower  than  t 

plasma  frequency.  The  impedance  probe  is  used  to  determine  values  o' 

electron  density,  electron  temperature,  and  collision  frequency  from  t n« 

difference  between  values  of  the  impedance  in  plasma  and  in  vacuum. 

Computer  developed  solutions  for  resistance  and  reactance  t v : ca 1 ! v  av- 

been  generated  and  graphed  as  a  function  of  plasma  frequency  and  t  <•!*>:•••?  , 

for  different  values  of  collision  frequency.  It  will  he  noted  '••a*  t  •:< 

primary  effect  of  electron  collisions  in  the  approximate  calculation  :>r< 

23 

cedures  developed  by  Despain  is  to  increase  the  antenna  resistance. 


The  data  gathered  by  the  Impedance  probe  allows  the  determination  o 
the  difference  of  the  antenna  impedance  from  free  space  impedance,  or 


AZ  «  AR  -  i  AX 

where  R  is  the  resistance  and  X  is  the  reactance  component  o'"  the  -.rwd 
Based  on  analysis  of  an  equivalent  circuit  for  the  probe,  the  follower 
equations  have  been  developed 

An  X  Z _ 

AR  =  o  'y 

to  c0[d  -  xr  +  ?-'j 


X(1  -  Xj _ _ 

u  Co((l  -  X)2  +  Z2] 


where  X  =  u)2  /to2  ,  i u j2  =  (n  e2/m  G  )  ,  Z  =  v/to  where  \  is  electron 
P  p  e  e  o 

collision  frequency,  C  =  C  +  C_  where  C  is  the  shunt  capacitance  or 

o  1  Z  o 
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ant  i-.'iru  miuDf  .iini  i  t  }  s  l  b»  tree  space  c.ifiar  i  t  ance  of  trie  antenna  . 
*•  t  r.t :  ehr  forward  fastium, 
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material  will  have  no  eft  '  on  n  or  .  m  • 1  •-  •>r<"  <  r.r  >• .  •• 

e 

the  uni  v  accurate  test  of  this  met  hod  won'.:  :nvol v  tl>«  f  .  •-  •  •  *  • 

and  so  determined  u  th  Independent  ,  e  era:  ••  n  .in'  v.t  1  ue  .  ••  •  .  •  ■ 

i 

c  omoa  r  :  son  o*  using  this  t ''c  hs.  i '■  ue  :s  no'  Known  to  t  c  , 

vein*,  can  he  nac**  that,  in  general,  n  values  derived  ‘  'or  ' «  f  <•  -:ri 

e 

are  higher  than  those  from  other  diagnostics.  The  ultimate  test  of 

method,  intercomparison,  does  seem  to  imnlv  a  difficulty  with  the  *  »-i  hn :  c  ut  . 

Differences  between  n  unleg  and  dovnleg  have  h'»"  ra  ’  ton. a  1  i  ?ed  T-  t  'ri  •  ■  a  ■* 

e 

q 

as  being  due  to  out.  gas ,  and  so  it  is  natural  to  examine  the  •ms 
influence  n>  out  gas  effects. 


First,  the  functional  forms  for  AR  (resistance)  and  AX  '  rea. 

demonstrate  opposing  effects  due  to  collision.  Consider  the  det  : na*  : r- 

of  X  from  AR,  C  ,  Z  or  AX,  C  ,  Z.  For  a  given  value  of  '.X,  inclusion  o' 
o  --  o 

outgas  enhanced  collision  effects  would  indicate  a  larger  np  than  without 

collision  effects.  For  a  given  value  of  AR,  however,  inclusion  of  outgas 

enhanced  collision  effects  would  indicate  a  smaller  n  than  without  collision 

-  e 


! 
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effects.  Recall,  however,  that  the  scheme  for  computation  of  does  not 

directly  include  a  collisional  model.  Without  a  new  fundamental  development 

of  impedance  relationships,  it  is  not  clear  whether  the  approximate  method 

of  n  determination  outlined  by  Despain  is  internally  self  consistent  with 

respect  to  collision  frequency.  Rather  than  to  judge  the  accuracy  of  that 

technique  with  respect  to  collision  freouencv  effects  on  n  determination, 

e 

nresent  efforts  will  concentrate  on  defining  physical  behavior  of  the  probe 
and  then  indicating  how  collision  effects  could  influence  the  determirj.it  ion 
of  n  .  gain,  collision  effects  produce  a  dominant  effect  on  the  real 
component  of  impedance,  the  resistance.  This  is  consistent  with  physical 
understanding  of  the  effects  of  electron  collisions,  i.e.,  there  is  a  direct 
reduction  of  effective  conductivity  of  a  medium  with  increased  ,  and  so  a 
dire 't  reduction  of  signal  amplitude.  Changes  in  phase  are  directlv  relat'd 
to  reactance  changes,  and  such  measurements  can  be  more  troublesome  to 
interpret.  Specifically,  if  one  takes  the  formula  for  X  as  a  function  o ‘ 

.  R  and  AX,  there  would  be  no  effect  on  n  due  to  outgas.  If  one  would  take 
the  formula  for  AX  =  f (X,  Z),  again,  there  is  little  effect  on  n^  determina¬ 
tion  due  to  Z.  However,  if  one  would  take  the  formula  for  .*  K  =  f  ?.'1  an<! 
use  it  to  determine  n  ,  there  would  be  definite  changes  in  'tv' i ra‘ ed  n  cv< 
*o  t  h*-  or'-since  of  outgas  material.  Some  details  will  *v>w  '->e  :  nee . 

The  effects  of  outgas  will  manifest  themselves  directlv  in  the  col  1  is 
frequency,  which  can  he  written 

1/2  r^Tp  1/2  r 

=  1.:'  i— -  1  ;  n.  Q  . 


IS 


where  Q  .  is  the  collision  cross  section.  For  the  altitudes  of  interest, 

e  j 

the  gases  are  primarily  1^*  so  we  can  write 


en 


{sec"1)  -  (f)1/2(^}1/2  {  (N2]Q 

e 

31  a,  ,  ,„-16  2 


2-eN2  +  [H2°lQeH20) 


Typical  values are  Q  „  %  4  x  10  10  cmZ  and  0  „  _  T  10  ^  cm*’  and  the 

eN2  eH.,0 

basis  for  anv  sensitivity  to  out  gas  can  be  seen  to  res-de  'n  the  coll is: or 
cross  section  which  :?  a  factor  of  10  higher  than  that  for  ir.  crcr 

to  gauge  effects  of  out  gas  in  this  fashion,  the  ca,f'.dat>n-  noted  above 
will  be  used:  a  constant  value  of  c  5  (  H,0 1  /  ’  N.,  ’  -  .  ?'.’4  v  '  ’  ’  V-  med  for 
all  altitudes.  In  order  to  estimate  the  electron  density  in  the  undisturbed 

medium  (without  outgas  or  shock  compression),  several  . .  must  he  taken . 

Values  of  electron  density  that  have  been  provided  were  deduced  by  straight¬ 
forward  application  of  X  =  f (AK,  AX).  However,  since  the  rocket  is  living 

sunerson icallv ,  correction  will  first  be  made  to  these  values  of  n  tc 

e 

account  *nr  both  shock  compression  and  electron  attachment  effects  (as  was 
done  for  the  supersonic  hangmuir  tip  probe)  to  attempt  to  recover  undisturbed 
electron  density.  These  will  be  presented  in  graphical  form  to  indicate  the 
magnitude  of  corrections.  A  second  procedure  to  account  for  outgas  effects 
will  then  be  carried  out.  It  initiates  by  computing  AR  from  X.  The 
effects  c‘  outgas  are  included  by  incorporating  Z  into  the  X  -  f (A?.,  7.) 
formula  to  get  a  new  (reduced)  n^  value.  Finally,  as  above,  the  effects  of 
shock  compression  and  electron  attachment  are  also  accounted  for  to  recover 
undisturbed  ,  in  this  second  case  with  corrections  for  outgas  and  shock 
effects . 

As  a  sample  calculation,  and  it  must  be  emphasized  that  data  from  a 

single  dav  is  limited  in  usefulness,  one  dav  on  which  comparisons  of  n  could 
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be  made  from  a  number  of  diagnostics  will  be  considered.  Data  were 
recorded  for  intercomparison  purposes  by  a  number  of  diagnostics  on  a  day 
preceding  the  eclipse,  but  part  of  the  1979  Canada-US  Solar 
Eclipse  Program.  Specifically,  data  were  recorded  by  a  blunt  probe  flown  by 
Mitchell  (UTEP),  supersonic  Langmuir  tip  probe  data  and  impedance  probe  data 
was  recorded  by  E.  Pound  (USU),  partial  reflection  data  were  recorded  by 
R.  0.  Olsen  (WSMR)  and  reduced  by  D.  Mott  (NMSU) .  The  presentation  of 

calculated  electron  densities  from  this  preliminary  data  is  first  made  in 
Figure  9.  The  blunt  probe,  tip  probe,  and  partial  indications  of  n0  have  been 
made  based  on  the  methods  noted  above.  The  Langmuir  tip  probe  data  was  cor¬ 
rected  for  shock  wave  and  attachment  effects  with  an  xl5  attachment  enhance¬ 
ment  factor.  The  impedance  probe  data  (ascent  and  descent)  is  presented  first 
without  any  adjustment;  second,  with  adjustment  for  shock  wave  compression  and 
attachment  effects;  third  for  outgas  adjustment  with  f  =  0.004  added  to  the 
shock  and  attachment  effects. 

The  presentations  of  calculated  electron  densities  including  all  the 
effects  and  techniques  noted  above  are  given  in  Figure  10.  For  this  partic¬ 
ular  day,  generally,  there  does  appear  to  be  reasonable  agreement  between 
supersonic  nose-tip  Langmuir  probe  indications  of  ng  and  those  from  the  imped¬ 
ance  probe.  Both  these  diagnostics  which  were  on  the  same  vehicle  do  differ 
above  65  km  from  the  values  deduced  from  blunt  probe  data  and  partial 
reflection  data.  It  is  noted  above  that  blunt  probe  analysis  above  65  km  is 
questionable  and  could  give  ne  values  that  are  up  to  an  order  of  magnitude  low 
at  75  km.  Also,  the  Inversion  of  the  partial  reflection  profile  above  65  km 
is  a  typical  output  when  the  technique  ceases  to  provide  accurate  ng  values. 
However,  these  were  very  active  days  and,  as  ng  point,  the  interval  of  time 
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Comparisons  of  defiled  calculations  of  electron  density  variations  for 
February  24,  19/9,  Red  Lake  (preliminary). 


Comparison  of  cal cul  ated  electron  density  variations  for  February  24,  1979, 
Red  lake. 


between  the  blunt  probe  data  shot  (1130  CST)  and  the  second  shot  with  the 
impedance  and  tip  probe  (1055  CST)  could  be  responsible  for  the  difference 
that  are  evident. 
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